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ABSTRACT 

We present simultaneous optical-to-X-ray spectral energy distributions (SEDs) from Swift's, 
X-ray and UV-optical telescopes (XRT and UVOT) for a well-selected sample of 26 low- 
redshift (z < 0.1) active galactic nuclei (AGN) from the Swift/Burst Alert Telescope (BAT) 
9-month catalogue, the largest well-studied, hard X-ray selected survey of local AGN to date. 
Our subsample consists of AGN with low intrinsic X-ray absorption (Ah < 10 22 cm~ 2 ) and 
minimal spectral complexity, to more accurately recover the intrinsic accretion luminosity in 
these sources. We perform a correction for host galaxy contamination in all available UVOT 
filter images to recover the intrinsic AGN emission, and estimate intrinsic dust extinction 
from the resultant nuclear SEDs. Black hole mass estimates are determined from the host- 
galaxy 2MASS K-band bulge luminosity. Accretion rates determined from our SEDs are on 
average low (Eddington ratios AEdd iS 0.1) and hard X-ray bolometric corrections cluster at 
~10-20, in contrast with the higher values seen for quasars. An average SED for the 22 low 
accretion rate (AEdd < 0.1) objects is presented, with and without correction for extinction. 
Significant dust reddening is found in some objects despite the selection of low Ah objects, 
emphasising the complex relationship between these two types of absorption. We do not find 
a correlation of optical-to-X-ray spectral index with Eddington ratio, regardless of the optical 
reference wavelength chosen for defining the spectral index. An anti-correlation of bolomet- 
ric correction with black hole mass may reinforce 'cosmic downsizing' scenarios, since the 
higher bolometric corrections at low mass would boost accretion rates in local, lower mass 
black holes. We also perform a basic analysis of the UVOT-derived host galaxy colours for 
our sample and find hosts cluster near the 'green valley' of the colour-magnitude diagram, 
but better quality images are needed for a more definitive analysis. The low accretion rates 
and bolometric corrections found for this representative low-redshift sample are of particular 
importance for studies of AGN accretion history. 
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1 INTRODUCTION 

Active Galactic Nuclei (AGN) are known to emit radiation over 
the whole range of available energies observable using current de- 
tectors. Characterising their spectral energy distribution (SED) is 
therefore important for understanding the different physical pro- 
cesses at work. The accepted paradigm behind the radiation output 
of AGN is accretion onto a supermassive black hole. The emission 
directly due to accretion emerges primarily in the optical, UV and 
X-ray regimes by a combination of thermal emission from an ac- 
cretion disc and inverse-Compton scattering of UV disc photons 
by a corona above the disc. A dusty torus partially absorbs the UV 
photons from the accretion disc, re-emitting them more isotropi- 
cally in the infrared. In this study, we aim to recover the true accre- 
tion luminosity in a representative sample of AGN, and use these 



results to identify trends between different SED parameters. How- 
ever, recovering the true accretion luminosity can be complex, and 
previous studies have found a wide spread in the conversion fac- 
tors (bolometric corrections) between the X-ray luminosity and the 
total accretion luminosity. This is, in a large part, due to intrinsic 
variation between objects in the whole AGN population, but is also 
complicated by numerous systematics and biases in the samples 
used. 

The pioneering study of lElvis et al.l d 1994) (hereafter E94) 
presented radio-to-X-ray SEDs for 47 quasars and has provided a 
useful template for determining bolometric corrections and iden- 
tifying trends in the SED shapes across a range of AGN prop- 
erties. Recent studies have confirmed the usefulness of the av- 
erage bolometric corrections and SED parameters from E94, but 
the sample of quasars used was predominantly X-ray bright and 
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their bolo metric luminos i ties i ncluded the re-processed infrared 
emission. iRichards et all d2006h present infrared-to-X-ray SEDs 
for 259 quasars selected by a combination of optical and mid- 
infrared colour selection criteria, and highlight the spread in SED 
parameters within the sample. This large spread in SED parame- 
ters was kno wn from E94, prompti ng refinements such as the tem- 
plate SED of iMarconi etafl d2004h who employ the observed cor- 
relation between X-ray-to-optical spectral index and X-ray lumi- 
nosity to construct a luminosity-dependent SED, and importantly, 
they exclude the reprocessed infrared emission to avoid double- 
counting of part of t he accretion energy b udget. A similar tem- 
plate is presented by iHopkins et ail d2007h in their study on the 
bolometric quasar luminosity function. These SED templates al- 
lowed the diversity of SED shapes to be taken into account when 
calculating the supermassive black hole mass density from the X- 
ray background and AGN luminosity functions. More recently, 
IVasudevan & Fabianl d2007l) (VF07 hereafter) presented optical- 
to-X-ray SEDs for a sample of AGN observed by the Far Ultra- 
violet Spectroscopic Explorer (FUSE), which yielded interesting 
trends between SED s hape a nd Eddington ratio, confirmed later by 
IVasudevan & Fabia n (2009) (VF09 hereafter) using simultaneous 
data and reverberation mapping mass estimates to improve accu- 
racy of the bolometric luminosities and accretion rates. However, 
the sample of VF07 is by necessity UV-bright, as it is selected 
by FUSE and the sample of VF09 is restricted to those which are 
bright enough in the optical/UV for reverberation mapping to be 
carried out. 

These considerations emphasise the pressing need for a rep- 
resentative sample of AGN from which conclusions can be drawn 
about accretion properties of the wider AGN population. One ma- 
jor factor influencing the selection of AGN for surveys is the nature 
of their absorption. It has been known for some time that X-ray and 
optical surv eys detect differen t segments of the underlying AGN 
population dMushotzkvll2004l) , with X-ray surveys able to probe 
to higher absorbing column density. The standard physical mech- 
anisms for the generation of X-ray emission in AGN produce a 
power-law, but often AGN exhibit X-ray spectra deviating signif- 
icantly from this, in part due to absorption. The geometry of the 
absorption in X-rays is the subject of much debate, and the com- 
plex spectra in some so urces can be accounted for with either 'par- 
tial covering' scenarios dGierliriski & D one 2004 ) or strong reflec- 
tion f rom the accretion disc due to light bending jMiniutti & Fabianl 
l2004h . and warm absorbers can also complicate our view o f the in- 
trinsic emission dRevnolds et al.lfl997l iBlustin et alj|2005h . Many 
of the reverberation mapped AGN in VF09 showed signs of such 
spectral complexity, and the authors identified the difficulties in re- 
gaining the true accretion luminosity in these cases. 

The Burst Alert Telescope (BAT) on board the Swift satel- 
lite proves invaluable for addressing many of these considera- 
tions. Canonical levels of absorption (column densities Nu < 
10 24 cm -2 ) imprint signatures on the 0.1-10 keV X-ray band, so 
the very hard X-ray sensitivity of the instrument (14-195 keV) pro- 
vides the capability to observe AGN in a bandpass relatively unaf- 
fected by such absorption. The 9-month catalogue of BAT-d etected 
AGN (hereafter the Swift/BAT catalogue. lfueller et alj2003) there- 
fore provides an unprecedented level of completeness (with respect 
to absorption) when surveying the AGN population as it is unbi- 
ased to all but the most heavily obscured sources. The study of 
IWinter et al] d2009r) presents a comprehensive overview of the X- 
ray spectral properties of the 153 sources in the Swift/BAT cata- 
logue and allows us to select an appropriate sample from which 
we can be confident of calculating accurate luminosities and SED 
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Table 1. Central frequencies of the UVOT niters. 

shapes. Most importantly, they present values of absorbing col- 
umn density Nn for each source, determined from fits to the X- 
ray data from the literature and new fits by the authors themselves, 
and identify those sources in which significant spectral complexity 
is present. In this work, we preferentially select those with lower 
absorption and minimal spectral complexity, facilitating a more 
straightforward calculation of the accretion luminosity. We assume 
that for such sources, the accretion luminosity excludes the repro- 
cessed IR emission and is principally seen in the optical-to-X-ray 
regime (as in VF07 and VF09). Additionally, AGN with more X- 
ray absorption are expected to display higher levels of optical-UV 
reddening, the calculation of which depends on the precise form of 
the extinction curve and is difficult to account for when calculating 
the total luminosity. These considerations motivate our selection of 
objects with low Nb_. 

Variability in AGN can produce large changes in luminosity, 
over timescales from hours to years. In order to catch an accurate 
snapshot of the total energy budget in an AGN at a given time, it 
is therefore highly desirable to use simultaneous data. The poten- 
tial inaccuracies in SED parameter values from non-simultaneous 
SED data are discussed in det ail in VF07 and VF09 , with the 
latter following the approach of iBrocksopp et all d2006l) by using 
contemporaneous optical-to-X-ray data from XMM-Newton's PN 
and Optical Monitor instruments. In this work, we make use of 
Swift's own co-aligned UV-optical telescope (UVOT) and X-ray 
telescope (XRT) to determine simultaneous SEDs. Many of our 
target subsample from the Swift/BAT catalogue have simultaneous 
UVOT and XRT observations available in the archives; typically 
the UVOT data span six filters between 5468 A and 1928 A, provid- 
ing a window onto a substantial fraction of the disc emission. The 
central frequencies of the UVOT filters are given in Table Q] 

The large wavelength coverage, field of view and spatial reso- 
lution of the UVOT data also allow a correction for the host galaxy 
flux to be estimated for each AGN. We employ GALFIT for this 
purpose, a program for performing 2D fitti ng of PSFs and va rious 
analytical profile forms to galaxy images dPeng etal.ll2002l) . For 
a detailed discussion of the requirements of GALFIT, the reader 
is directed to peruse the accompanying documentatiorQ. We also 
present basic results on the host galaxies from the GALFIT fitting 
process. The UVOT data used in this study also clarify how the de- 
gree of host galaxy contamination varies with wavelength (UVOT 
filter). 

Understanding variations in the AGN SED is important, and 
one of the main drivers of these variations appears to be the ac- 
cretion rate (parameterised as the Eddington ratio, Lboi/iEdd 
for an object with bolometric luminosity Lboi and Eddington lu- 
minosity L E dd = 1.38 x 10 38 [MBH/Mo]ergs _1 for a black 
hole of mass Mbh). The recent Principal Component Analysis of 

1 http://users.ociw.edu/peng/work/galfit/galfit.html 
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iKuraszkiewicz et al.l < l2008t) on a sample of red 2MASS AGN SEDs 
identify the Eddington ratio as the 'eigenvector 1' for these AGN 
(being responsible for most of the variation in the SED shapes). 
Both VF09 and VF07 highlight substantial variations in the bolo- 
metric output (in comparison to the X-ray output) across the ob- 
served range of Edd ington ratios, and numerous studies ex plore 
these themes further iShemmer etai]|2008l. [Kelly et al.ll2008l) . Ed- 
dington ratios require accurate estimates of the black hole mass 
in AGN, but since the most accurate estimates from reverbera- 
tion mapping are only available for 35 AGN, here we employ the 
Mbh — -t-buige correlation to calculate black hole masses. The com- 
bination of robust determinations of the bolometric luminosity and 
sensible estimates of the black hole mass provide great opportu- 
nities for insight into the accretion process. In this work, we dis- 
cuss the analysis of the Swift data, the construction of SEDs, the 
derivation of SED parameters, the identification of correlations and 
trends, and finally the host galaxy properties. 



2 SAMPLE SELECTION 

We firstly apply a cut in intrinsic column density N^ nt \ select- 
ing only those with log(iv4 lnt ')< 22 from the list of ob j ects i n 
the 9-month BAT catalogue as presented in IWinter et alj | |2009|) . 
This is the suggested crossover point between 'absorbed' and 'un- 
absorbed' classes of AGN in their study, which also classifies the 
objects in the sample into two broad categories based on X-ray 
spectral properties: 'simple' (S) and 'complex' (C). The S class 
have X-ray spectra which are best fit by a simple power law with 
intrinsic and galactic absorption (sometimes with a soft excess), 
whereas the C class display more complex X-ray spectra which are 
better fit by models such as double-power laws or partial cover- 
ing (these models are illustrative in that they highlight the spectral 
complexity; more detailed fits including reflection may also be pos- 
sible). We also eliminate any objects from the C class, as they will 
present difficulties when determining the true accretion luminosity 
(as discussed in VF09). 

In order to facilitate profile fits to optical and UV images of 
the host galaxies, we also impose a redshift cut, requiring z < 0.1 
for our selection. For z < 0.1, the angular resolution of the UVOT 
(~ 1.7 arcsec) allows separation of the nucleus from the galaxy 
to a physical length scale of 3 kpc, thus allowing a reasonable 
galaxy-AGN separation while still yielding 54 potential objects 
for study from the sample. We then identify those for which Swift 
XRT and UVOT data are available from the High Energy Astro- 
physics Science Archive Research Center (HEASAItd^), yielding 
33 objects. For estimating black hole masses, we also obtain K- 
band magnitudes from the Two-micron All-Sky Survey (2MASS) 
catalogues. The key requirement for calculating black hole masses 
is an estimate of the bulge luminosity. The K-band is least subject 
to the effects of Galactic reddening and predominantly traces the 
older stars in the bulge over the stellar populations of the galaxy 
disc, motivating its selection over the J and H bands also avail- 
able in the 2MASS catalogues. The 2MASS data were gathered 
with two ground based telescopes: the Whipple Observatory in Ari- 
zona, USA and the Cerro Tololo telescope at La Serena, Chile. The 
limitations introduced by ground-based observations (most impor- 
tantly, the level of seeing for each observation) need to be taken 
into account when attempting to recover the bulge luminosity. This 

2 http://heasarc.gsfc.nasa.gov/cgi-bin/W3Browse/swift.pl 



is discussed at greater length in i]3.4l on black hole mass determi- 
nation. We required the centroids reported for the 2MASS data to 
be within 2.0 arcsec of the NED position in the initial case. For 
a few objects which lie just outside this limit, the 2MASS images 
were downloaded and the locations provided were checked against 
the refined positions obtained from the XRT and UVOT images; if 
a confident identification of a counterpart to the X-ray AGN was 
made, the 2MASS data were used in further processing. At any 
rate, objects with 2MASS centroids greater than 5.0 arcsec away 
from their NED positions were rejected altogether. Three potential 
objects are excluded from our sample by their absence from the cat- 
alogues (IGR 21277+5656, IGR J07597-3842 and LEDA 138501). 

Lastly, we check for blazars in the sample. The SED proper- 
ties of blazars make them unsuitable for performing a bolometric 
luminosity calculation based on the 'disk + power-law' model em- 
ployed here since their UV and X-ray luminosity is dominated by 
jet emission. Out of the 30 objects which satisfy all the above cri- 
teria, four objects are exlcuded based on their NED classifaction as 
blazars: Mrk 501, 2MASX J19595975+6508547, ESO 362-G021 
(PKS 0521-36) and PKS 0548-322. The final selection of 26 ob- 
jects is presented in Table[2] along with their redshifts and details 
of the Swift observations used. 



3 PROCESSING OF MULTIBAND DATA 

We download the optical-UV and X-ray data for our selection from 
HEASARC. When multiple observations were available, those 
with the maximum UVOT exposure time and maximum number 
of UVOT filters were preferentially selected. Pipeline-processed 
'level 2' FITS files are readily available from HEASARC. The data 
for the 26 sources identified in section [2] were then processed ac- 
cording to the procedure outlined in the following sections. 

3.1 Optical-UV photometry: Correcting for the host galaxy 

We employ the custom-built software tools designed specifically 
for processing UVOT data where possible. Each individual UVOT 
filter data file in general contains a number of exposures which 
were summed using the tool UVOTIMSUM. The UVOTSOURCE tool 
was then used to extract magnitudes from simple aperture pho- 
tometry. Source and background regions were created for this pur- 
pose, with the position of the source region being obtained from 
the NASA Extragalactic Database (NED) in the first instance, fol- 
lowed by fine adjustment of the source region position if necessary. 
The required source region size for UVOTSOURCE is 5 arcsec, and 
the background region can be any size (at the time of writing): we 
used background regions with radii from 5 arcsec to ~30 arcsec de- 
pending on the frequency of other foreground sources in the image. 
The magnitudes from UVOTSOURCE provide a useful first-order 
estimate of the nuclear flux, and were saved for a later compari- 
son with the more carefully determined nuclear fluxes using point 
spread function (PSF) fitting. 

The images were then prepared for use within GALFIT. One 
of the requirements for accurately determining the nuclear flux is 
to have suitable PSFs available for each image. Preliminary stud- 
ies of the UVOT PSF in the six filters show that the PSF shape 
may depend on a variety of factors including position on the de- 
tector, countrate of the source and the filter being used. Average 
PSF full-width half-maxima range between 1.7-2.5 arcsec, w hen 
considering all six filters. As discussed by iKim et all d2008ah in 
the context of Hubble Space Telescope (HST) images, the choice 
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Table 2. Details of UVOT observations used for the objects in our sample. 



of PSF can significantly affect PSF fitting results, since a poorly 
chosen PSF can give profile para meters that deviate dram atically 
from their 'true' values (also see ISimmons &Urrvll2008l for an- 
other discussion of the factors affecting AGN-host-galaxy decom- 
position). To mitigate these problems and minimize PSF mismatch, 
we adopted the following approaches. The first and preferred ap- 
proach was to generate a unique PSF for each filter image from a 
few known guide stars in the image, requiring the stars selected to 
be within a certain range of the count rate of the source of inter- 
est (as determined from simple aperture photometry with UVOT- 
SOURCE). Firstly, the UVOT utility UVOTDETECT was used to pro- 
vide a list of detectable sources in each image (UVOTSOURCE es- 
sentially uses the SEXTRACTOR package for this purpose). These 
were then cross-checked against published catalogues of known 
guide star positions. The XBROWSE tool in the HEASOFT suite of 
utilities provides access to many such catalogues, such as the US 
Naval Observatory (USNO) and HST guide star catalogues. Typi- 
cally, the USNO catalogue was used to obtain a list of such stars, 
requiring them to be within 13 arcmin of the AGN of interest to 
ensure they were within the UVOT field of view. The guide stars 
were selected such that between 3-20 stars were used to generate 
each PSF image. For the lowest energy V-band, this translated into 
a requirement that the stars were typically within 0.1-0.2 dex of 
the target AGN count rate; for the highest energy UVW2 band, the 
limits were typically within 1.5-2.0 dex of the target AGN count 
rate to obtain a similar number of guide stars. Finally, the regions 
of the images containing these stars were then summed to form 
a final PSF image, using IDL code by A. A. Breevelcjf] designed 
specifically for generating PSFs from UVOT images. The resulting 



PSFs were used in GALFIT to model the central AGN in the galaxy 
profiles. 

When this approach was not possible due to lack of guide 
star detections within reasonable count-rate limits (or lack of point 
sources which were also identified in the guide star catalogue), the 
filter image was viewed using the XIMAGE package and a 'King' 
profile model was fit to an identifiable point source which looked 
similar in form and intensity to the AGN nucleus. The parameters 
from the fit were then used as a model for the PSF in GALFIT (via 
the 'Moffat' model option), in lieu of a real PSF image. In the vast 
majority of cases, a real PSF image could be extracted from the im- 
ages and a 'King/Moffat' model PSF was rarely required. In only 
one filter image for one object, a Gaussian model was found to fit 
the nucleus better than other available models, with the Gaussian 
width fixed at the average FWHM of the UVOT PSF in that filter 
as reported in the CALDB documentatiorjfl 

The PSF or PSF model thus generated was then fit along with 
other galaxy profile components in GALFIT. We performed tests 
to ascertain the level of detail discernible within the UVOT im- 
ages, starting initially with a four-component model consisting of 
a constant sky background, a central nuclear point source, a Ser- 
sic profile for the galaxy bulge and an outer exponential disk. Such 
tests overwhelmingly indicated that four components incorporated 
too high a degree of degeneracy, due to the large UVOT PSF size. 
The PSF of the UVOT instrument is typically the same size as the 
galaxy bulge would be in higher resolution images of galaxies near 
redshift ~0.05, and indeed for some of the more distant sources, 
the UVOT PSF is as large as the host galaxy itself. Often the Ser- 
sic profile component would shrink to a size comparable to that of 



3 Mullard Space Science Laboratory - http://www.mssl.ucl.ac.uk/ 
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the PSF, and would 'absorb' some of the flux which would other- 
wise been reported as coming from the PSF component. Therefore, 
three components were adopted as sufficient to give sensible model 
reconstructions of the profiles seen in the images (sky + PSF + ex- 
ponential disk). To get as accurate a representation as possible of 
the nuclear flux, the following iterative algorthim was employed 
using GALFIT: 

(i) All bright foreground stars were excluded from the fit by cre- 
ating a bad pixel mask which covered these unwanted objects. 

(ii) The sky background was fit independently to a blank region 
of the sky close to the AGN of interest. 

(iii) The central point source and sky were fit together, keep- 
ing the sky parameters constant from the previous step. The central 
point source magnitude was seeded with that obtained from UVOT- 
SOURCE aperture photometry, to provide a sensible initial guess of 
the nuclear flux. 

(iv) The sky, central source and exponential disk components 
were fitted together, keeping the sky and central source parameters 
fixed at the values from the previous step. 

After these steps, the resultant model fits were compared with 
the original images and the residual images obtained from subtract- 
ing the model from the data were inspected. In the first instance, 
the positions of the different components were left free to fit. If any 
obvious signs of PSF mismatch or other problems with the fit were 
evident, these positions were frozen manually by inspection and 
the fit was re-evaluated. In some cases, such as for point-like AGN 
at the higher end of our redshift range, a simpler model consisting 
of just a PSF and a constant background was adopted, if the first- 
pass attempt at fitting an exponential disk did not show a significant 
galaxy disk component. Where the PSF image displayed obvious 
PSF mismatch, a 'King' profile fit was also attempted. The best fit 
obtained from all these approaches, determined by visual inspec- 
tion of the residual image, was chosen to obtain the final nuclear 
magnitudes. An example of the sky-nucleus-disk decomposition is 
shown in Fig. [2] The choices of profiles used for each filter in each 
object are given in table[3] 

3.1.1 Comparison with HST 

The detailed study of lBentz et al. highlights the need to take 

into account host galaxy contamination when calculating nuclear 
luminosities from optical HST images, specifically in the context 
of determining the radius-luminosity relationship for the broad line 
region (BLR) in AGN. They find that even for an aperture of 1 arc- 
sec, there can be significant host galaxy contamination. However, 
the HST Advanced Camera for Surveys (ACS) offers a significant 
advantage with a PSF of FWHM 0.0575 arcsec in contrast to the 
typical FWHM of ~2 arcsec for the UVOT PSFs. This somewhat 
limits our ability to determine the host galaxy contamination, since 
the stars used to model the PSF often are as wide in extent on the 
sky as some of the more distant galaxy disks in our sample. We 
attempt to gauge the possible underestimate of host galaxy contri- 
bution with UVOT using publically available images from the HST 
ACS and Wide Field and Planetary Camera 2 (WFPC2) for one of 
our objects, NGC 4593. This object is known to have a clearly iden- 
tifiable spiral galaxy structure surrounding the central AGN so pro- 
vides a good foundation for comparison between UVOT and HST 
image quality. Our UVOT analysis also shows a clear discrepancy 
between the 5.0 arcsec aperture photometry results and those from 
PSF fitting in the optical (V, B and U) bands. Since some level of 
galaxy contamination is clearly identifiable using UVOT images, 



it should therefore be more readily identifiable in the HST images. 
We perform simple aperture photometry using the DS9 package on 
images from ACS for the F550M and F330W filters (centered on 
5580A and 3353A respectively). We extract counts from the nu- 
cleus with a circular region with radius 0.2 arcsec, which should 
provide an upper limit on the nuclear flux and contain >90 per cent 
of the counts from the nucleus. We also determine the total counts 
from a 5 arcsec region, as used with UVOTSOURCE for comparison. 

In the F550M filter (closest to the UVOT V-band in wave- 
length), we identify a difference in fluxes of a factor ~5 between 
the smaller and larger apertures for HST. In comparison, the differ- 
ence in fluxes between the 5.0 arcsec UVOTSOURCE magnitude and 
that from PSF fitting to the UVOT image is a factor ~2.4. Assum- 
ing the 5.0 arcsec fluxes from the two instruments should theoreti- 
cally be identical, this implies that the flux reported to the nucleus 
by PSF fitting to UVOT images is a factor ~ 5.0/2.4 too large. The 
problem is less acute in the F330W filter (closest to the UVOT U- 
band), with an identical analysis yielding that the nuclear flux from 
the UVOT image is a factor ~ 1.8/1.2 of that from the HST im- 
age. There is also the issue of optical AGN variability, which could 
be artificially increasing these ratios (the continuum near 5100A 
is known to vary by factors of ~ 7 in the case of NGC 5548; see 
IPeterson et al.ll 19991) , but it is difficult to disentangle this from the 
error intrinsic to the instrument, and the factors we calculate here 
are therefore probably upper limits. This reflects a fundamental 
limitation of the resolution of the UVOT images, but the subse- 
quent analysis shows that in many of the objects seen, attempting 
to remove the host galaxy does produce an optical-UV SED shape 
closer to that expected for an AGN despite possible overestimates 
of the nuclear flux. A more detailed analysis using observations 
from the Kitt Peak National Observatory (KPNO) 2.1m telescope 
should provide a much more detailed galaxy profile decomposition 
(Koss et al, in prep), which will make use of the smaller PSF avail- 
able with the KPNO (FWHM ~ 1.0-1. 4 arcsec). 

3.2 Optical-UV photometry: Correcting for coincidence loss 
and Galactic extinction 

For AGN with very high count rates, the fluxes obtained from 
UVOT images are highly susceptible to coincidence losses and 
therefore often appear lower than their true values. Coincidence 
loss refers to the phenomenon where multiple photons arrive at the 
same location on the detector durin g a single frame, and is dis- 
cussed in detail in lPoole et al.l j2008l) . They the theoretical correc- 
tion to be applied to the measured count rate for individual pixels, 
but also point out the need to account for the spread of a point 
source over many pixels. They also provide an empirical correc- 
tion which takes the latter effect into account. We correct all of 
our magnitudes from GALFIT for coincidence loss manually using 
these expressions; the UVOTSOURCE aperture photometry magni- 
tudes are automatically corrected for coincidence loss. 

There is some scope for error in the coincidence loss cor- 
rection when the PSF image is not quite a point source, as the 
coincidence loss correction only holds good for point sources. 
This problem is especially pronounced for very bright sources. In 
these sources, the coincidence-corrected fluxes are sometimes even 
greater than those obtained from UVOTSOURCE, which should the- 
oretically provide an upper limiting nuclear flux. In cases where 
this phenomenon is very pronounced, we do not use the magni- 
tudes from GALFIT and instead employ the 5.0 arcsec aperture 
photometry magnitudes in any further analysis. 

After calculation of magnitudes from both GALFIT and 
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Figure 1. Example GALFIT profile fitting results for MCG-06-30-15. The first column in each row shows the original data image, the second shows the model 
fit and the third shows the residual obtained from subtracting the model from the data. The rows, from top to bottom, show the results for the V, B, U, UVW1, 
UVM2 and UVW2 bands, in that order. In the V and B band, the residuals imply that a bulge component may also be discernible in this object, but the UVOT 
data did not allow robust fits including a bulge component for many of our objects. A simple analysis using DS9 regions indicate that the putative bulge may 
account for up to 10 per cent of the total galaxy flux in the V and B band images. 



UVOTSOURCE, these magnitudes were corrected for Galactic ex- 
tinction using the values for E(B — l / )cai from the NASA/IPAC 
Infrared Science Ar chivjfl and the Galactic extinction curve of 
ICardelli et af] ( ll989h . The final magnitudes were turned into XSPEC 
PHA files using the FLX2XSP utility, to facilitate fitting along with 
the contemporaneous XRT data using the XSPEC analysis software. 



http://irsa.ipac.caltech.edu/applications/DUST/ 



3.3 X-ray spectra from XRT 



The pipeline-processed event files from the XRT detector were pro- 
cessed using the XSELECT package, as directed in the Swift XRT 
user guide. Source regions of 50 arcsec were used, with larger 
accompanying background regions (average radius ~ 150 arcsec). 
Background light curves were determined from the event files and 
inspected for flaring, but this was not found to be a problem in any 
of our observations. Source and background spectra were extracted, 
and the source spectra were grouped with a minimum of 20 counts 
per bin. 
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AGN 



UVWl UVM2 UVW2 



1RXS J045205.00+493248 


P 


P 


P 


P 


_ 




2MASX J21 140128+8204483 


P+E 


P 


P+E 


P 


P+E 


P+E 


3C 120 


P+E 


P+E 


P+E 


P+E 


M+E 


P+E 


3C 390.3 


P+E 


P+E 


P 


P 


P+E 


P+E 


Ark 120 


P+E 


P+E 


P+E 


P+E 


P+E 


P+E 


ESO 490-G026 


P+E 


P+E 


P+E 


P+E 


P+E 


P+E 


ESO 548-G081* 


P+E 


P+E 


P+E 


P+E 


M+E 


M+E 


IRAS 05589+2828 


P 


P 


P 


P 


M 


Gaussian 


IRAS 09149-6206 


P+E 


M 


P+E 


P+E 


P+E 


P+E 


MCG +04-22-042 


P+E 


P+E 


P+E 


P+E 


M+E 


P+E 


MCG-06-30-15 


P+E 


P+E 


P+E 


P+E 


P+E 


P+E 


Mrk 1018 


P+E 


P+E 


P+E 


P+E 


M+E 


M+E 


Mrk 279 


P+E 


P+E 


P+E 


P+E 


P+E 


P+E 


Mrk 352 


P+E 


P+E 


P+E 


P+E 


P+E 


M+E 


Mrk 509 


P+E 


P+E 


P+E 


P+E 


P+E 




Mrk 590 


P+E 


P+E 


P+E 


P+E 


P+E 


P+E 


Mrk 766 


P+E 


P+E 


P+E 


P+E 


P+E 


P+E 


Mrk 841 


P+E 


P+E 


P+E 


P+E 


M+E 


P+E 


NGC 4593 


P+E 


P+E 


P+E 


P+E 


P+E 


P+E 


NGC 5548 


P+E 


P+E 


P+E 


P+E 


P+E 


P+E 


NGC 7469 


P+E 


P+E 


P+E 


P+E 


P+E 


P+E 


NGC 985 


P+E 


P+E 


P+E 


P+E 


P+E 


P+E 


SBS 1136+594 


P 


P 


P+E 


P+E 


P+E 


P+E 


SBS 1301+540 


P+E 


P 


P+E 


P 


P 


P+E 


UGC 06728 


P+E 


P+E 


P+E 


P+E 




P+E 


WKK 1263 


P+E 


P+E 


P 


P 


P 


P 



Table 3. Table showing the model components used in GALFIT with each filter image for each AGN. 'P' signifies a PSF generated from guide stars within 
appropriate count-rate ranges, from each image, 'M' signifies a Moffat (King) profile with parameters estimated from a single point source in the image, and 
'E' signifies an exponential disk model component.. *For ESO 548-G081, an extremely bright foreground star was located within ~30 arcsec of the AGN, 
rendering the photometry obtained susceptible to large errors. 





1000 2000 3000 



Figure 2. Comparison between images from HST (left) and UVOT (right) for the source NGC 4593. A 5 arcsec circular region is shown for reference on both 
images. 



3.4 Black hole mass estimates using the 2MASS All-Sky 
Survey Catalogues 

The Mbh — ibuigc correlation for galaxies provides a useful way 
of estimating their central black hole masses. The key challenge is 
to obtain an accu r ate est imate of the bulge luminosity. The study of 
iMarconi & Huntl d2003l) presents a decomposition of images from 
2MASS using GALFIT to obtain the bulge luminosity and corre- 
late it with the black hole mass from direct determination meth- 
ods (stellar or gas kinematics, maser kinematics, etc.). However, 
the redshifts of the objects u nder scrutiny here are ge nerally sig- 
nificantly higher than those in lMarconi & Hunl d2003h . presenting 
problems for per f ormin g a full nucleus-bulge-disc decomposition. 
iMushotzkv et al.l d2008h present a method for calculation of black 



hole mass for the Swift-BAT catalogue AGN from 2MASS K-band 
total source magnitudes, which involves subtracting the central nu- 
clear component from the extended source flux to estimate the 
bulge luminosity. We attempt to refine this method here by incorpo- 
rating information on the expected angular size of the bulge on the 
sky and the resolution limitations of the 2MASS data. The crucial 
factor is the degree of seeing affecting the 2MASS observations; 
typically the seeing takes a value of 2.6 arcsec . For b ulges with 
typical size ~ 0.5 - 1 kpc dGraham & Worley||2008l) . the bulge 
becomes unresolvable at redshifts of about z ~ 0.01 (assuming 
Ho = 71 kms -1 Mpc" 1 , = 0.27 and a flat Universe). There- 
fore, for the overwhelming majority of objects in our sample, the 
bulge should be unresolved. Inspection of the atlas images from 
2MASS for the lowest redshift objects in our sample confirms that 
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the visually identifiable bulge component is not resolved. The bulge 
light is therefore mixed in with the nuclear light in these sources 
and identified as a 'point source' . The point sources identified by 
2MASS have been collated, along with their magnitudes, in the 
2MASS Point Source Catalogue (PSC); similarly extended sources 
have been catalogued in the 2MASS Extended Source Catalogue 
(XSC). Based on our assumption that the bulge is unresolved, we 
initially download only the magnitudes from the PSC for further 
use, as they will contain the bulk of the bulge light we wish to re- 
cover. 

By way of calibrating our attempt to obtain black hole masses, 
we first turn to the reverberation m apping (RM) sample of AGN 
presented by IPeterson et alj d2004l) . representing the most secure 
mass deter minations for AGN t o date (augmented with the refined 
result from lDennev et al1l2006l for NGC 4593). There are uncer- 
tainties in the RM method, connected to assumptions about the ge- 
ometry of the broad-line region and whether it is gravitationally 
bound, but the ~ 35 AGN for which RM mass estimates exist nev- 
ertheless constitute the sample of AGN with the most carefully de- 
termined black hole masses. Initially, we download the 2MASS at- 
las images for those RM AGN below a redshift of 0.01, in which 
the bulge should be theoretically resolved and perform a GALFIT 
3-component decomposition using a gaussian PSF for the nuclear 
point source, a Sersic bulge and an exponential disc. Our GALFIT 
analysis yields that there is a significant degree of uncertainty in 
the identification of the bulge component in these AGN, and the 
magnitude ascribed to the bulge is heavily dependent on the prior 
constraints imposed. We explore a number of different constraints 
on the model profiles; for example, the Sersic index n is often con- 
strained to an upper limit to avoid the fitting routine confusing it 
with the nuclear point source, and often requires a lower limit to 
avoid confusion with the disc component. Despite such strategies, 
in some sources the bulge and disc profiles acquire very similar or 
identical radii and 'share' the disc luminosity, with the Sersic pro- 
file partially tracing the disc instead of the bulge. In such sources, 
three-component model fits do not yeild meaningful bulge parame- 
ters. Our attempts at constraining the bulge using GALFIT result in 
bulge luminosities that vary by up to an order of magnitude depend- 
ing on the choice of model and constraints imposed, and show that 
even for the close by sources, the bulges are probably not resolved. 

In light of this, we then employ the 2MASS PSC catalogue 
magnitudes to provide a simple estimate of the total bulge and nu- 
clear flux. In order to estimate the relative contributions of the nu- 
cleus and bulge in the PS C magnitudes, we employ the infrared 
SED templates presented in lSilva et alj d2004h . In their study, they 
present nuclear IR SED templates constructed using available near- 
to-mid IR data on 33 Seyferts, using a radiative transfer models for 
dust heating to interpolate between the wavelengths covered by the 
data. They also present host galaxy SED templates by subtracting 
the nuclear template contribution from the total photometry at each 
available wavelength. They present a number of different nuclear 
SED templates appropriate for different levels of X-ray absorption, 
and host SED templates are presented for different intrinsic 2-10 
keV AGN luminosity (Lx) regimes. Using the appropriate nuclear 
SED along with the corresponding host galaxy SED, it is there- 
fore possible to estimate the fraction of the total luminosity which 
can be accounted for by the host for an AGN with a particular Lx- 
In the process we make the assumption that, in the case of the K- 
band, the host flux should be dominated by the bulge as discussed 
in iJT] We calculate the ratio of host to total luminosity (nuclear 
plus host luminosity) from these SED templates in the K-band for 
the four different luminosity regimes (41.0 < log(Lx) < 42.0, 



42.0 < log(Lx) < 43.0, 43.0 < log(Lx) < 44.0, 44.0 < 
log(Lx) < 45.0). The fractions are plotted in Fig. [3] against the 
central luminosity of each bin. We also plot a logarithmic inter- 
polation between the points, to allow determination of the ratio 
of Lk, host / Lk, total at a general value of Lx, assuming the ratio 
varies continuous l y. For values of Lx greater than those spanned 
by the lSilva et all 020041) study, we simply employ the fraction cal- 
culated for Lx = 10 45 ergs _1 , extrapolating from the direction 
of the trend shown between the bins centred on log(Lx) = 43.5 
and 44.5. For luminosities log(Lx) < 41.5, we use the fraction 
evaluated at log(Lx) = 41.5 (extrapolating from the trend would 
yield fractions larger than Unity). The high luminosity extrapola- 
tion is clearly not likely to be accurate for extremely powerful ob- 
jects (log Lx > 46), but none of the reverberation mapped AGN 
or the AGN in our Swift-BAT subsample are in this regime. In the 
absence of information on IR SEDs outside this luminosity range 
we do not attempt a more complex extrapolation. 

Before using these fractions to calculate bulge luminosities 
from the 2M ASS PSC magnitud es, we consider potential sources 
of bias. Since [silva et alj | |2004|) adopt a reprocessing scenario to 
account for the K-band continuum, it is possible that the model nu- 
clear K-band luminosity is a function of the Eddington ratio (since 
it is ultimately linked to the UV accretion disc luminosity, the dom- 
inant component of the total accretion luminosity). The X-ray lu- 
minosity can be recast as the product of the Eddington ratio, bolo- 
metric correction and black hole mass (via the Eddington ratio), 
so the trend seen in Fig. [3] could be the result of two functions 
of Eddington ratio being plotted against each other. However, we 
again note the low X-ray luminosities probed by the Swift/BAT 
9-month catalogue (up to a few times 10 44 ergs _1 ), and find that 
the host-to-total fraction varies at most by a factor of two in this 
regime. Therefore, the dependence is not strong and will not intro- 
duce significant biases compared to other sources of error, such as 
the intrinsic spread in RM masses and the intrinsic dispersion of 
the Mbh — Lk, bulge relation. Regardless of the particular model 
used to predict the nuclear K-band continuum, it is clear that the 
fraction of an unresolved bulge that is attributable to the nucleus 
will decrease for lower luminosity AGN and vice versa; the simple 
relationship in Fig. [3] provides an observationally-rooted estimate 
of this effect. 

We employ the v alues of Lx from XMM-Newton reported in 
IVasudevan & Fabianl j2009t) to calculate the K-band host-to-total 
ratio for each of the RM objects, and assume the I Silva et alj d2004l> 
Seyfert 1 IR SED template for all objects in the sample. In the ab- 
sence of a detailed, uniform study on the absorbing columns of the 
RM AGN sample this is a reasonable assumption, since the rever- 
beration mapping sample is heavily biased towards Seyfert 1 AGN; 
none of the RM AGN are classified as Seyfert 2s in NED. We scale 
the PSC magnitude by the appropriate fraction, using the expres- 
sion 

M K ,bulgc = M K ,PSC - 2.51oglo (/bulge (ix)), (1) 

where the fraction /bui ge = Lk, host /Lk, total is calculated as 
above. The black hol e masses are then calcula ted using the Mbh — 
Lk, bulge relation of iMarconi & Huntl J2003h . The resulting mass 
estimates are compared with RM mass estimates in Fig. [4] 

A reasonable agreement between the two methods is seen, al- 
beit with some significant deviation at higher masses. The grey 
shaded area highlights the uncertainties in RM masses due to the 
lack of precise knowledge of the BLR geometry. The largest dis- 
crepancy is for the AGN with the highest mass, 3C 273, for which 
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Figure 3. Variation of Lj<, host /-^K, total with X-ray luminosity, extrap- 
olated from the host and nuclear IR SED templates of ISilva etal.1 {2004) 
(black filled points - determined using the four fiducial luminosities pro- 
vided in their paper; solid line - extrapolation between points as detailed 
in the text). These fractions are employed to calculate the contribution of 
the bulge to the 2MASS point source magnitudes, for determining black 
hole masses. The curve for an AGN with intrinsic X-ray absorption Nn = 
10 23,5 cm~ 2 is shown (dashed line) for comparison (the curves for other 
absorptions above Nn = 10 22 cm — 2 are very similar; the significant dif- 
ference is between obscured and unobscured AGN). 



the K-band bulge luminosity estimate gives a black hole mass an 
order of magnitude too high. This is known to be a very power- 
ful object, and it is possible that our simple approach to estimating 
the bulge fraction in the point source light has still underestimated 
the powerful nucleus. There is also likely to be a significant syn- 
chrotron (non-accretion) contribution to fluxes at all wavelengths 
in this source, which further increases uncertainties. This AGN is 
far more powerful than the majority of the Swift-BAT catalogue 
AGN, so such problems are not likely to affect the lower-power 
obj e cts we consider here. 

iMarconi et al. point out that the previous reverberation 

mass estimates have been calculated with the assumption that the 
BLR is gravitationally bound, and introduce a correction based on 
an estimator of the accretion rate (the 5 100 A monochromatic lu- 
minosity) to provide revised mass estimates for the RM AGN cat- 
alogue. We also provide a comparison with their revised mass es- 
timates, in Fig. [5] The correlation between the two mass estima- 
tors is strengthened significantly if these first-order corrections to 
the reverberation masses are applied, and particularly reduces the 
scatter seen for those objects with high accretion rates (for which 
the correction is maximal and works to increase the reverberation 
mass estimate). The offset is still present however, and a simple fit 
of the form log(M B H,L K bulgc ) = A + log(M B H,rcvam P ) yields 
A ~ 0.57 in both cases (using reverberation masses with or with- 
out correction for radiation pressure), implying an offset of a factor 
of up to ~ 3.7 in black hole mass, just outside the typical quoted 
tolerance for reverberation masses (a factor of ~ 3). We now sug- 



Figure 4. Comparison between masses from reverberation mapping and 
masses estimated from K-band bulge luminosities. The solid line represents 
the desired one-to-one correspondence between the two methods and the 
shaded area depicts the intrinsic uncertainties in reverberation mapping. The 
dashed line shows the best fit of the form log(MeH,L K bu igo) = A + 

log(M B H, revamp)- 



gest a possible contributor to this offset and present a correction for 
it. 

The lSilva etal] J2004l) host galaxy SEDs have been calcu- 
lated using large-aperture and extended-source data (such as the 
2MASS XSC catalogue magnitudes) which could include more 
tha n the bulge light. I f the galaxy disc is significant in some cases, 
the lSilva et al.l d2004l ) host- to-total ratios actually provide the frac- 
tion of bulge plus galaxy disc to the total luminosity rather than 
just the bulge, and therefore using them to scale the 'nucleus plus 
bulge' luminosities from the 2MASS PSC is not strictly correct. It 
is therefore necessary to gauge the effect of neglecting the galaxy 
disc in this calculation. We do this by obtaining the 2MASS Ex- 
tended Source Catalogue magnitudes (matching the centroids ap- 
propriately as before), using these fluxes as estimates of the total 
nuclear, bulge and disc light (X t otai = L nuc + L bulgc + L disc ). 
Employing the PSC magnitudes as the bulge and nuclear light as 
before (Lb u i gc + -Lnuc), we then can ut ilise the fra c tion o f nuclear- 
to-total luminosity predicted from the ISilva et all d2004l) IR SED 
templates (/ nuc = L nuc /[L nuc + L bulgc + L disc \) to obtain a the- 
oretically more accurate estimate of the bulge magnitude: 



M K.PSC _M K,XSC 
2.51ogl (l - /nuclO ■ O )), 

(2) 

This expression is obtained by solving for Lbuigc from the avail- 
able quantities. Put simply, it arises from using the nuclear frac- 
tion from the IR SED templates to estimate the nuclear flux from 
the XSC flux; this nuclear component is then subtracted from 
the PSC flux to obtain the bulge flux. This estimate is not avail- 
able for all of the objects available previously, as some objects 
do not have 2MASS XSC magnitudes. For a couple of sources, 
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M BH from reverberation mapping (corrected for radiation pressure effects) 



Figure 5. Comparison between masses from RM ( corrected for radiation 
pressure effects as discussed by Marconi et al. 2008) and masses estimated 
from K-band bulge luminosities. The empty white circles show the results 
using the reverberation masses before correcting for radiation pressure. 
Other key conventions as in Fig. [4] 

the pred icted fraction / nu c = L nuc / (L nuc + Lb u i ge ISC / 

from the lSilvaetafl fe004h templates was larger than the fraction 
(Lnuc+-Lbuigc)/(£nuc+Lbuigc+£di S c) found from the ratio of the 
PSC and XSC fluxes, leading to an impossible solution for Lbuige- 
This indicates tha t the nuclear-to-total fraction obtained from the 
ISilva et alj J2004h templates is not appropriate for that object, or 
that the bulge is not completely contained in the PSC flux. These 
problems limit the use of this method, but for those objects where a 
calculation is possible, we find that the new K-band mass estimates 
are systematically shifted down by 0.14 dex (a factor 0.72) from 
the values from the simpler approach of using the PSC magnitude 
alone. We therefore apply this correction to all masses obtained 
using the PSC alone, which brings the best-fit line for the new 
mass estimates within the tolerance of the reverberation masses. 
The strength of the correlation is good, as indicated by the Pear- 
son correlation coefficients of 0.83 and 0.93 obtained (correlating 
with reverberation masses before and after radiation pressure cor- 
rection). This analysis also indicates that the galaxy disc contribu- 
tion in the host-to-total ratios is of the order ~ 30 per cent, and a 
similar value applies to all AGN in the RM sample. There still re- 
mains some offset above reverberation and within the tolerance re- 
gion, but any further attempt to account for this is somewhat redun- 
dant in light of the intrinsic uncertainties in reverberation masses, 
so we do not make any other corrections to our approach. The off- 
set of virial mass estimates below the Mbh — Lbuige r elation has 
been noted and discussed in detail bv lKimetalJj2008bh : they find 
that it is a complex function of a variety of parameters including 
the black hole accretion rate and host galaxy morphology. 

We finally calculate the masses using this approach for all of 
the z < 0.1 objects in the Swift-BAT 9-month catalogue for which 
2MASS PSC magnitudes are available. The K-band host-to-total 
ratios are calculated using the Lx values reported in lWinter et al.1 
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Figure 6. Comparison between masses from reverberation mapping and 
masses estimated from K-band bulge luminosities (with offset due to un- 
derestimation of the galaxy disc contribution taken into account). Key as in 
Fig.|5] 

d2009h . and the A^h values in the same paper are used to select be- 
tween the two host-to-total ratio curves for obscured or unobscured 
AGN (Fig. [3} as required. By way of another check on the accu- 
racy of our mass estimates, the recent study o f IWang et al.l ( l2009h 
presents virial mass estimates for a sample of Seyferts which par- 
tially overlaps with Swift-BAT 9-month catalogue AGN, using the 
single-e poch H/3 linewidth- based black hole mass estimator (as 
given in iGreene & Ho1l2005l) . We plot a comparison between their 
values and ours in Fig.|7]for the 24 overlapping objects. While the 
relation shows some scatter, we caution that the linewidth-based 
black hole mass estimation methods also show significant scatter 
themselves when compared to RM. The two estimators still appear 
to trace each other (with a correlation coefficient of 0.80); addi- 
tionally, preliminary calculations of black hole masses from Kin 
Peak National Observatory (KPNO) and Sloan Digital Sky Survey 
(SDSS) spectra for the Swift/BAT 9-month catalogue AGN also 
seem to match our estimates well, with a scatter of a factor of ~ 2 
(Winter et al., Koss et al. in prep). We see some indications that 
the local Swift-BAT 9-month catalogue AGN have predominantly 
low black hole masses, as expected in 'anti-hierarchical' structure 
formation scenarios. The masses thus determined are used in all 
subsequent analyses for calculation of accretion rates and for con- 
straining accretion disc models in SED construction. 



4 GENERATION OF SPECTRAL ENERGY 
DISTRIBUTIONS 

We fit the optical, UV and X-ray SED data points with a simple 
multicolour disc DISKPN and absorbed power-law model combi- 
nation using XSPEC, following the approach outlined in VF09. In 
this work, we extend the approach of VF09 by inclusion of a multi- 
plicative ZDUST component on the DISKPN model used for the UV 
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Figur e 7. Comparison between masses from H — j3 linewidths iWang et"aT] 
120091) and masses estimated from K-band bulge luminosities. The solid line 
repr esents the d e sired one-to-one correspondence between the two meth- 
ods. [\Vang|etal] b009l) do not provide error estimates on their black hole 
masses, and the typical (random) error in the K-band estimates can be 
gauged from Figs [4] and [5] (it is typically dominated by uncertainties in- 
trinsic to the Mbh — ^bulgc relation), although the discussion in the text 
highlights that systematic errors are likely to dominate. 

bump, to estimate intrinsic reddening where possible. The presence 
of data in at least four UVOT filters for all of the objects in the sam- 
ple makes this possible, and indeed since most of the objects have 
UVOT data in all six filters, we are able to constrain the shape of 
the SED over a wider energy range than was generally possible in 
VF09. Signs of intrinsic reddening (such as an optical-UV spec- 
trum steeply declining with energy) can then be more unambigu- 
ously identified. 

The full model used in XSPEC consisted of the optical- 
UV reddening component multiplied by a thermal accretion disk 
model, and a broken power law in X-rays with both Galac- 
tic and intrinsic absorption components (ZDUST(DISKPN) + 
WABS(ZWABS(BKNPOWER))). The normalisation of the DISKPN 
component was determined as in VF09, assuming an inclination 
angle of zero and a colo ur-to- temperature ratio of Unity. Unifica- 
tion scenarios for AGN l lAntonuccill 1993b suggest that AGN with 
lower gas column densities such as the ones in this sample should 
be preferentially oriented with their accretion disks almost face-on, 
and therefore small inclination angles are appropriate (see VF09 for 
a discussion of the effect of varying inclination angle in the DISKPN 
component on luminosities calculated from SEDs). We also freeze 
the inner radius of the disc at 6.0 gravitational radii, which is appro- 
priate for radiatively efficient accretion. The local Galactic absorp- 
tion was determined using the NH utility from the FTOOLS suite 
and frozen in the fit. The low energy branch of the BKNPOWER 
model was set to have negligible values within the energy range 
of the UV bump (as in VF09). We fit the XRT data in the range 
0.3-10.0 keV. While there is an interesting debate as to the pre - 
cise nature of the intrinsic reddening in AGN jGaskell et ai]|2004l) . 



we adopt the Small Magellanic Cloud (SMC) reddening curve for 
our 'first-order' estimation of intrinsic extinction E(B — V) and 
freeze the ratio of total to selective extinction, Ry, at 3.0. The in- 
trinsic extinction E(B — V), the intrinsic column density N^ nt \ 
the maximum temperature T max for the DISKPN model and the pa- 
rameters for the high-energy branch of the BKNPOWER component 
were all left free in the fit. Occasionally, if the combined optical, 
UV and X-ray data resulted in a photon index in the X-ray regime 
which was influenced by the slope between optical and X-ray, the 
photon index was frozen at the value obtained from just fitting the 
X-ray data alone. 

We freeze the normalisation of the UV bump DISKPN compo- 
nent using the black hole mass estimate from the 2MASS K-band 
bulge luminosity estimate (normalisation K = Af| H cos(i) / Z?l/3 4 
for black hole mass A/bh, inclination angle i, luminosity distance 
Dl and colour-to-effective temperature ratio (3). We present the 
SEDs determined using this approach in Fig. [8] For a few objects, 
the DISKPN model thus constrained fits poorly to the data despite 
their relatively 'blue' optical-UV spectral shapes, thought to be 
typical for the lower-energy part of an un-reddened accretion disc 
spectrum (i.e. Ark 120, Mrk 1018, Mrk 509 and NGC 985), which 
is probably attributable to a slightly too large black hole mass. This 
is confirmed later in the case of Ark 120 for which a reverbera- 
tion mass estimate is available (see §53\ . In the case of Mrk 590, 
the distinctly reddened optical-UV spectral shape is poorly fit us- 
ing the 2MASS black hole mass-derived DISKPN normalisation, but 
the reverberation mass estimate produces a more plausible fit. The 
DISKPN model fit for IRAS 09149-6206 also looks poor, although 
this object possesses a rather un usual, flat optical-UV SED (noted 
initially bv lKirhakos et al.ll 199 lh unlike the other AGN in the sam- 
ple, which may be due to some other process. These five cases 
might point to a genuine offset between the reverberation masses 
and our 2MASS derived estimates, but the degeneracy of the mass 
with the inclination and other parameters in the DISKPN normal- 
isation could indicate, for example, unusual inclination properties 
in these few AGN. The remaining 21 objects all exhibit sensible 
model fits to the optical-UV points. 

We present the average model SED (before and after correc- 
tion for reddening) calculated for the 22 AGN with Asdci < 0.1 
(normalised at leV) in Fig. [9] This allows a useful comparison with 
the average SEDs presented in VF07 and VF09 for high and low 
Eddington ratio AGN, and the mean quasar SED templates of E94; 
the latter are also reproduced in Fig. [9] While the E94 template is 
representative of X-ray bright quasars, ours offers a similar tem- 
plate for local, X-ray unabsorbed and low accretion rate Seyferts, 
albeit calculated from the model fits instead of the data. Some dif- 
ferences between our average SED and those of E94 include the 
slightly reddened shape of the E94 SEDs in the optical-UV (their 
SEDs trace our extincted SED more closely, at least at longer wave- 
lengths) and the differing X-ray spectral shapes. The radio-quiet 
quasar X-ray spectrum in particular appears softer in the 0.5-2 keV 
band which would be expected for a brighter, high er accretion rate 
samp le of AGN (according to the trends seen in e.g. lShemmer et all 
|2008|); but this could also be in part due to the inclusion of some 
X-ray absorbed objects in the E94 composite, depressing the flux at 
this part of the spectrum relative to the flux at higher energies. The 
hard shape of the radio-loud template can be accounted for by the 
expected synchrotron contributions from jets. The E94 SEDs also 
show a prominent soft excess; due to the limited spectral resolution 
of the XRT we do not consider the soft excesses in our SEDs here. 
We repeat the cautionary note from E94, that AGN properties de- 
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Figure 8. A selection of AGN SEDs. The crosses represent the 5 arcsec aperture photometry results from UVOTSOURCE, and the filled circles represent the 
data points used for the fit. In the UV and optical regime, the black filled points represent the fluxes with the host galaxy component removed using GALFIT 
as detailed in the text. Empty circles were not used in the fit in a few cases where significant coincidence losses risk affecting the accuracy of the fluxes. 



termined from the average SED do not necessarily reflect the full 
degree of variation in the AGN population. 

We follow the approach outlined in VF09 for calculating the 
unabsorbed hard X-ray (2-10 keV) and bolometric (0.001-100 
keV) luminosities and errors in XSPEC, for use in determining ac- 
cretion rates, bolometric corrections and other SED parameters. 
The results are presented in Table IaTI As discussed in fJT] we de- 
fine the bolometric luminosity as being the primary emission due 
to accretion emerging in the optical-to-X-ray regime and assume 
that the IR is reprocessed emission; it is therefore not included in 

Lbol- 



4.1 Extending the SEDs into hard X-rays 

The BAT data allows us to extend the SEDs into hard X-rays (14— 
195 keV). Four-channel BAT spectra are publicly available on the 
interne^, time-averaged over periods of many months. We iden- 
tify the four least variable objects in the 14-195 keV band using 

B http://swift.gsfc.nasa.gov/docs/swift/results/bs9mon/ 



measurements of excess variance on the BAT lightcurves, in order 
to maximize confidence in luminosities extrapolated from fits to 
the combined X-ray the hard X-ray data. We present the SEDs for 
these four sources (NGC 5548, NGC 7469, Mrk 279 and 2MASX 
J21 140128+8204483) including the 4-channel BAT data in Fig.fTOl 
We also present the fits to the more robust 8-channel data (the data 
themselves are not yet public at the time of writing) to illustrate 
any substantive differences between 4- and 8-channel datasets. In 
all cases, the simple model fits (POWERLAW or PEXRAV) do not in- 
dicate a signifcant departure from the general shape implied by the 
4-channel data. 

In three out of four of these low-variability sources, it ap- 
pears that the BAT data may be capturing a reflection hump 
peaking at around ~30 keV, in accordance with standard sce- 
narios for reflection from the accretion disc. There does not 
seem to be any pronounced reflection signature in NGC 5548. 
We corroborate this suggestion by fitting the model combination 
[WABS(ZWABSPOWERLAW)+PEXRAV] to the XRT and BAT data 
together in XSPEC, linking the photon index of the illuminating 
source in PEXRAV to that for the POWERLAW component. The 
model fits (thin dot-dashed lines) agree well with the fits to the 
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Figure 10. Spectral energy distributions with BAT data included. Empty squares represent 4-channel BAT data (publicly available). Thick dashed lines 
represent fits to the 8-channel data (not public at the time of writing). The thin dot-dashed line represents a fit to the XRT and BAT data combined, using the 
model [WABS(ZWABSPOWERLAW)4-PEXRAV] to illustrate the plausibility of a reflection fit to the data. 
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Figure 9. Average, reddening-corrected optical-to-X-ray SED for the 22 
low Eddington ratio (Agcid < 0-1) AGN, normalised at 1 eV (thick black 
solid line). The average fit (in the optical-UV regime) with reddening in- 
cluded is shown using the thick black dashed line . The mean SEDs for 
radio-quiet and radio-loud quasars from lElvis et alj 1 19941) are also repro- 
duced (thin black solid and dot-dashed lines, respectively). 



XRT alone at low energies, and show a plausible reflection hump 
fit to the BAT data at higher energies for NGC 7469, Mrk 279 and 
2MASX J21 140128+8204483. There are hints of accompanying Fe 
K-a lines at 6.4 keV in the XRT data, as expected from reflection 



models, but the quality of the XRT data does not lend itself to a 
confident identification of this line. The superior quality XMM-PN 
data presented for NGC 7469, Mrk 279 and NGC 5548 in VF09 
show a much more convincing Fe K-a line shape at 6.4 keV, al- 
though the line signature is less evident in NGC 5548, which would 
be expected for weaker reflection. The indications of reflection in 
these sources are confirmed by a more detailed analysis using XMM 
and BAT data (Mushotzky et al. 2009 in prep). 



5 DISCUSSION 

5.1 Contamination of optical-UV continuum with emission 
lines 

We briefly comment on the effect of emission lines in optical-UV 
spectra of AGN which have the potential to int roduce systematic ef- 
fects in our SED measurements. The study o f lFrancis et alU 19911) 
presents a high signal-to-noise composite quasar spectrum, allow- 
ing the parameters of such emission lines to be accurately deter- 
mined. We identify the emission lines with the largest widths and 
highest peaks from Figs. 2 and 7 of their work and redshift them in- 
dividually for each AGN in our sample to ascertain whether emis- 
sion line fluxes are likely to contaminate the UVOT photometry. 

As a first-order attempt to quantify this, we ascertain whether 
the central wavelength of the redshifted line lies within the band- 
pass of each of the UVOT filters, as obtained from the CALDB 
documentation. The most common contaminant from this anal- 
ysis would seem to be the Mg II line at 2798A (rest-frame), 
which is most likely to appear in the UVW1 filter, and oc- 
casionally in the U filter for two high-redshift objects. Al III 
(195 8 A) and C III] (1909 A) are also likely to be common con- 
taminants in UVM2 and UVW2 bands. For the four highest red- 
shift sources (3C 390.3, IRAS 09149-6206, SBS 1136+594 and 
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2MASX J21 140128+8204483), the redshifted O III lines (4959A 
and 5007 A) appear in the UVOT V-band and in the case of 2MASX 
J21 140128+8204483, the H/3 line (4861 A) also. The presence of 
these lines may explain, in part, the 'angular' behaviour of the 
optical-UV SEDs in these four sources. However, the V-band 
'bumps' in more local sources (ESO 548-G081, Mrk 279, NGC 
4593 and NGC 7469 for example) are not accounted for by the 
presence of such lines and must have some other cause. In any 
case, for the vast majority of objects in the sample, the SED points 
in most filters are highly unlikely to be affected by the presence of 
any prominent emission lines. 

5.2 Intrinsic reddening due to dust extinction 

We present the values for intinsic extinction E(B — V) determined 
from the SEDs against Eddington ratio in Fig. QT] This allows a 
comparison between the amount of estimated dust extinction in the 
source and the accretion power available from the central AGN. 
One notable feature of the distribution of E(B — V) observed is 
that despite the selection of objects with low Nn, there appears 
to be significant dust reddening in many of these sources, based 
on the sha pe of their optical-UV spectrum. This was noted pre- 
viously by iKraemer et al. I d2000h and provided the basis for their 
'lukewarm absorber' scenari o to account for this di screpancy. Our 
results contrast with those o f iMaiolino etail ( 1200 ll) . who find that 
dust-to-gas ratios - E(B — V)/Nh - are systematically lower in 
their sample AGN than in the Galaxy. However, their sample is 
selected specifically to include higher Njj. objects with relatively 
un-absorbed optical/UV/IR broad lines, and the data used are not 
simultaneous. The results of both of these studies results imply that 
objects which are significantly absorbed in X-rays can show lit- 
tle sign of dust extinction in the optical-UV and vice-versa, and 
in general that classifications of AGN into 'obscured' or 'unob- 
scured' based on their optical properties can differ significantly 
from X-ray classifications. Our sample of objects confirms that op- 
tical and X-ray obscuration scenarios can differ significantly in an 
object, influen ced as they are by different physical processes (as 
summarised in ICrenshaw & KraemejlioOlh . We caution of a de- 
generacy between the mass and E(B — V) in the model fits (not- 
ing in the case of Mrk 590 that using its RM mass gives a best-fit 
E(B - V) of 0.33 compared to ~ 0.1 using the 2MASS mass 
estimate), although the SEDs for ESO 490-G026, MCG-06-30-15, 
Mrk 590 and Mrk 766 show more unambiguous signs of significant 
intrinsic extinction. In support of this, the aforementioned stu dy on 
reddened Seyfert 1 galaxies of ICrenshaw & Kraemed d200lh con- 
firms that observations from the International Ultraviolet Explorer 
(IUE) also strongly suggest intrinsic reddening in MCG-06-30-15, 
Mrk 590 and Mrk 766. The high o ptical polarisation for Mrk 766 
discussed bv lUlvestad et al.l i 19951) may also be an indicator of the 
dust responsible for reddening in this sour ce. Encouragingly, the 
more detailed study of MCG-06-30-15 by iRevnolds et all Il997h 
reports a lower limit on E(B — V) very similar to the value of ~0.6 
that we determine for this source, implying that our estimation of 
intrinsic extinction with the ZDUST model can provide reddening 

estimates in line wit h more detailed stud ies. 

The studies of iFabian etafl d2008h and iFabian etail J2009h 
highlight the connection between the radiation pressure exerted by 
the AGN and the configuration of the surrounding absorbing ma- 
terial. They identify the effective Eddington limit for dusty gas 
in the Nu — Asdd plane; above this limit is a 'forbidden region' 
in Nn — AEdd space within which absorbing dusty gas clouds 
are unstable to radiation. They verify that most local and deep 



0.6 



0.5 



9 



I 0.3 
.= 

a 0.2 



0.1 





1 1 

— •MCG -06-30-015 


-^Mrk 766 




5 09149-6206 

O 490-G026 , 

»BGC 06728 

^Mrk 352 


t*nggXL 


-WKK 1263 Isrs H3 6+5 94 
-f3C 390.3T_^ SBS | 301+540 

f Ml 



10" 3 



0.01 0.1 
Eddington ratio 



Figure 11. Intrinsic extinction E(B — V) against Eddington ratio Asdd> 
for those objects for which reddening was included in the model fit via the 
ZDUST XSPEC model. 



samples of AGN avoid the forbidden region, and the latter work 
in particular verifies this in the case of the unbiased Swift/BAT 
catalogue. Objects near the effective Eddington limit are known 
to exhibit AGN winds/warm absorbers, which would be expected 
when radiation pressure locally exceeds gravity. Since the dust 
and gas are expected to be coupled, the dust extinction parame- 
terised by E(B — V) may also exhibit a similar distribution in 
the E(B — V) — Asdd plane. We see from Fig. II II that the ma- 
jority of objects occupy the lower-left part of the plot, and the 
only object with particularly h i gh red dening - MCG-06-30-15 - is 



i particularly higl 

also known from I Winter et al .1 ( 120090 to have a relatively high A"h 
(~ 1.9 x 10 21 cm -2 ) and AEdd, placing it close to (but not within) 
the forbidden region. This object showcases a warm absorber, as 
expected for an object near the effective Eddington limit. The ab- 
sorbing column reported from analyses of Beppo SAX and ASCA 
data i s strongly model-dependent (see for example. IRevnolds et al .1 
1 19971 who find Nu to be an order of ma gnitude lower than the Win- 
ter et al. result; or lMorales et alj|2000L whose detailed multi-zone 
warm absorber scenario gives a column density of 3 x 10 22 cm~ 2 ). 

These hints may suggest that an analogous 'forbidden region' 
could be proposed in the upper-right portion of the E(B — V) — 
AEdd plane. A more complete study of intrinsic reddening proper- 
ties for a larger sample of sources spanning a wider range of X-ray 
column densities would be valuable for exploring this possibility 
further. 



5.3 Bolometric corrections and Eddington ratios 

We present the bolometric corrections K2-iokcV = 
iboi/i2-iokcv against Eddington ratio in Fig. [T2] The bolo- 
metric corrections cluster between 10-20 with a low fraction of 
objects possessing higher bolometric corrections, and Eddington 
ratios are overwhelmingly below < 0.1. The low bolometric 
corrections obtained are expected for low Eddington ratios (as 
found in VF07 and VF09) and the few objects with high bolometric 
corrections generally lie at higher Eddington ratios as expected. 
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Errors in luminosities take into account the errors in model 
parameters obtained from fitting the data, and the Eddington ratio 
additionally takes into account the random error in the black hole 
mass estimate (extrapolated from the uncertainties in the K-band 
magnitudes). However, the potential systematic errors due to 
intrinsic disparities between RM and K-band mass determination 
methods also contribute an uncertainty on both axes which is not 
represented in the error bars. 

We therefore invoke RM ma ss estimates for nine s ources com- 
mon to the objects discussed bv lPeterson et alj d2004l) to perform 
some simple comparisons which provide an estimate of the magni- 
tude of this systematic effect. If the RM masses are used in the fit- 
ting process instead, we obtain bolometric corrections and Edding- 
ton ratios as given in Fig.[T3] We also present the results from VF09 
for comparison, which employed simultaneous data from XMM- 
Newton. If reddening is not included, as in VF09, we obtain very 
similar bolometric corrections and Eddington ratios as seen with 
XMM (left panel, Fig. 1131). However if reddening is taken into ac- 
count, this generally increases both bolometric corrections and Ed- 
dington ratios (right panel, Fig. 1 1 3t and in the case of objects such 
as Mrk 590 and NGC 7469, this change is very pronounced. How- 
ever, in both scenarios, the general trend for bolometric correction 
increasing with Eddington ratio identified in VF07 and VF09 is 
preserved. 



7469 _ 



4 Mrk5 "° - f" ^ , 2. 

y -^Ark 120 

-^C 120 

10-1 ^«^? C5MS >C3M3 

590 -4mGC 5548 . 

- T [ »3€ 390.3 



0.01 0.1 1 

^Edd 

Figure 14. Hard X-ray (2-10keV) bolometric correction against Edding- 
ton ratio for the objects with reverberation mapping mass estimates. Here 
we present a comparison between the values obtained using mass estimates 
from 2MASS K-band bulge luminosities (blue empty squares) and reverber- 
ation mapping (black), using the UVOT and XRT data in both cases, with 
intrinsic reddening included in the fit. 



We also perform a comparison between RM masses and 
masses from K-band luminosities using only the Swift simultaneous 
data, including intrinsic reddening effects in both cases (Fig. 1 14b. 
Bolometric corrections and Eddington ratios are clustered at signif- 
icantly higher values for this subsample when using the RM masses 
compared to the K-band masses. There is an increase of around ~ 1 
dex in Eddington ratios when RM masses are used, accompanied 
by an increase in bolometric correction to an average of ~ 30. The 
most extreme changes are for NGC 7469, Mrk 590 and Mrk 279. 
These results illustrate the degree of change which could be seen if 
the 2MASS K-band mass estimates need to be scaled down to bring 
them more in line with RM estimates, but given the BLR geome- 
try uncertainty intrinsic to the RM masses, the bolometric correc- 
tions and Eddington ratios calculated using better-calibrated RM 
masses could also be reduced significantly. As depicted in Fig. [6] 
the two methods are broadly consistent within these uncertainties. 
We suggest that the true magnitude of any shift in results may lie 
in between these two extremes; the 2MASS masses may represent 
an upper limit whereas the reverberation estimates may represent a 
lower limit. Indeed, the radiation pressure corrections m ay require 
further calibration as discussed bv lMarconi et aljj2009l) . and could 
bring these estimates closer together, for example. More work is 
required on calibration between different black hole mass estima- 
tors for AGN before this issue can be addressed more convincingly, 
and at this stage we are only able to provide a crude estimate of the 
limiting uncertainties. 

These analyses suggest that the typical bolometric corrections 
and Eddington ratios for this sample are low, even when uncer- 
tainties in the mass estimates are taken into account. The distribu- 
tion of Eddington ratios and bolometric corrections are presented in 
Fig.[T5] using the 2MASS K-band mass estimates. Th e distribution 
is sim ilar to that presented for Seyfert Is in Fig. 7 of Iwinter et al.1 
j2009t) shifted by approximately an order of magnitude, in line with 
the relatively low bolometric corrections we find for the majority of 
the sample. 



5.4 Fraction of ionizing luminosity 

The ionizing continuum in AGN is responsible for the production 
of the optical and UV emission lines which are often an integral 
part of selection criteria in optical surveys of AGN. Variations in 
the ionizing fraction can therefore lead to pronounced changes in 
the detection of these lines (VF07). We plot the ionizing fraction 
(£e>t3.6ev/£boi) against X-ray luminosity, bolometric luminos- 
ity and Eddington ratio in Fig.[T6]to identify any possible trends. 
It appears that there is little correlation between Le>iz.6bV / L\> \ 
and luminosity (either X-ray or bolometric), but there may be indi- 
cation of a more substantial correlation with Eddington ratio. This 
is in line with the variation in SED shape reported in VF07 and 
VF09, despite our sample spanning only the lower-Eddington ratio 
regime. The average ionizing fraction is ~0.6, consi stent with that 
adopte d for the ionising luminosity corretion b lon in lMarconi et al.1 
d2008l) . but there is a large spread; for low Eddington ratio objects 
this could drop to ~0.4, and rise to ~0.8 for high Eddington ra- 
tios. We note that these ionizing fractions are extrapolated from the 
DISKPN big blue bump model fits, and that there may be significant 
components contributing to this ionizing flux in the extreme-UV 
and soft X-ray regime, which our approach cannot quantify; the 
ionizing luminosity fractions presented here may therefore repre- 
sent lower limits to the 'true' ionizing fractions. However, such a 
trend for increasing ionizing fraction with Eddington ratio is also 
expected from combined ac cretion disk and corona models such as 
the one bv lWitt et al.Nl997h . The distrubution of values seen for the 
ionizing luminosity fraction may also need to be taken into account 
in the radiation pressure corrections to virial masses as discussed 
bv lMarconi et al.1 J2009h . in a refinement over their 2008 work. 

5.5 Optical-X-ray connections 

The Qox-Liy ( 2500A ) correlation report ed in the literature 
dStrateva et all I2005L ISteffen et al.1 l2006h has been widely 
used to encapsulate a relationship between the UV accre- 
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Figure 12. Hard X-ray (2-10keV) bolometric correction against Eddington ratio. The mass estimates used for model fitting and calculation of Eddington ratios 
are determined from the 2MASS bulge luminosity. 
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Figure 13. Hard X-ray (2-10keV) bolometric correction against Eddington ratio for the objects with reverberation mapping mass estimates (black filled 
circles). Left panel: Intrinsic reddening was not included in the fit, as was the case in VF09. Right panel: Intrinsic reddening was included in the fit, to contrast 
with VF09 and assess the role that intrinsic extinction may play in shaping the SED for some of the reverberation mapped AGN. In both cases, the reverberation 
mapping masses were used in the fit and to determine the Eddington ratios. The results from XMM-Newton PN and Optical Monitor (VF09) are shown for 
comparison (blue empty squares). 
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Figure 15. Histograms of Eddington ratios (left panel) and bolometric correction (right panel) for the sample. 



42 42.5 



43 43.5 
log(L 2 „| 0keV ) 



44.5 45 45.5 

iogtu.0 




Figure 16. Ionizing fraction £_E>l3.6ev/^bol- Left panel: L E>i3.6eV I '-t<bol against X-ray luminosity. Middle panel: iB>i3.6ev/^bol against bolometric 
luminosity. Right panel: L E>i3.6eV / Lbol against Eddington ratio. 



tion disk emission and the coronal X-ray emission, and the 
relation between optical and X-ray emission in AGN has 
been e xtensively studied observationally (e.g. the RIXOS cata- 
logue, iMason et alj l200d) . The spectral index aox is defined 
as -log(L„(2500A)/L^(2keV))/log(i/(2500A)/z/(2keV)) and 
therefore the optical flux used in its calculation is measured at 
an energy significantly lower than the peak of the big blue bump 
(in the vicinity of ~ 1000 A for typical AGN black hole masses of 
~10 7 ~ 8 Mq, accreting at ~10 per cent of the Eddington limit). We 
assess how changing the optical reference point in the calculation 
of this optical-to-X-ray spectral index alters the correlations seen, 
and denote the new spectral indices as oibx , otux and oiuvwix 
(which use the rest-frame B, U and UVW2 bands as the optical 
reference point instead of 2500 A). It is evident from the center- 
left panel of Fig. ( 1 7) that this sam ple follows the established cor- 
relation from lSteffen et alj | |2006|) . within the expected scatter. We 
present spectral indices calculated both from the raw data (by in- 
terpolating to the desired optical wavelength, black filled circles) 
and from the de-reddened model fit (white filled circles). In the 
case of aj3x, there seems to be little correlation as the B-band is at 



the low-energy end of the disk component, where variations in the 
disk emission have relatively little effect. However the correlations 
linking qux and Quvw2x to their respective reference luminosi- 
ties seem more substantial for the raw, interpolated values: the cor- 
relation coefficients decrease monotonically from —0.46 to —0.79 
as the reference wavelength moves from the B to the UVW2 band, 
demonstrating the stronger anti-correlations seen using a shorter 
reference wavelengths. However, the intrinsic spectral indices do 
not show such trends: the correlation coefficients actually increase 
from —0.27 to —0.16 going from B to UVW2, which on visual in- 
spection of these plots reveals a large degree of scatter between the 
reddening-corrected spectral indices and their reference luminosi- 
ties. 

This demonstrates the importance of accounting for redden- 
ing. Previous studies on aox have attempted to minimize the ef- 
fects o f reddening by carefu lly selecting the samples used; for ex- 
ample IStrateva et alj J2005t) require that their AGN are required 
to be bluer than a template AGN SED with fixed reddening; red- 
shifting this template then provid es a redshift-dependent colour 
threshold (see lRichards et alj2003l) which they use to exclude 'red' 
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Figure 17. From left to right, top to bottom: Spectral indices connecting SED at the rest frame B-band, U-band, 2500A («ox) an d UVW2-band with the SED 
at rest frame 2keV. Each spectral index is plotted against the corresponding rest-frame reference luminosity in the optical/UV band (e.g. Lb for obx> L\j 
for ctuxetc). The black filled circles represent the values obtained from interpolating at rest-frame energies from the data points (not corrected for intrinsic 
reddening) and the white filled circles s how the indices calc ulated from the de-reddened SED model fit. The solid line and shaded area in the top right panel 
shows the best fit and spread obtained bv lSteffen et alj bfJOd) for a sample of 333 optically selected AGN. The correlation coefficients of de-reddened, intrinsic 



spectral index against luminosity are as follows: obx — Lb>0.27, aux — Ly. -0.24, aox ~ L 



-0.19, «UVW2X — Luvw2>0.16. 



AGN. This should statistically exclude heavily reddened AGN, 
but we find from individually de-reddening our optical-UV AGN 
SEDs, that reddening effects could still be significant. The anti- 
correlations between the raw spectral indicies (uncorrected for red- 
dening) and their respective reference luminosities appear strong 
because reddening in an object produces a flatter optical-X-ray 
spectral index and reduces the UV luminosity; this will tend to pop- 
ulate the top-left regions of each plot with reddened objects while 
the lower-right regions will contain unreddened objects, produc- 
ing an artifically robust-looking relationship. This effect may con- 
tribute to the strength of the aox — ^ 2 500A correlat i° n m tne liter- 
ature. From the perspective of understanding the accretion process, 
however, it is the intrinsic, de-reddened optical-X-ray spectral in- 
dices which are of interest. We only use a small sample here, but 
if individually de-reddened optical SEDs for AGN can be used to 
construct the aox — L ■ relation anew with larger samples, it 
will clarify whether the versions of this relation seen in the litera- 
ture can be used to robustly constrain the accretion process in AGN. 

When plotted against Eddington ratio instead of optical/UV 
reference luminosity, even the anti-correlations between the raw, 
uncorrected indices and optical/UV reference luminosities disap- 
pear (see Fig. 1 1 81>. This phenomenon was also observed by VF07 



and lShemmer et all J2008h . The latter study suggests that the lack 
of significant correlation between optical-X-ray spectral index and 
Eddington ratio could be because of inaccuracies in mass determi- 
nation methods. They also discuss how the use of an X-ray selected 
sample (such as the hard X-ray selected BAT catalogue used here) 
should alleviate some selection effects which artificially could blur 
any correlation with Eddington ratio, but this study seems to rein- 
force the lack of connection between optical-X-ray spectral index 
and Eddington ratio. This seems to be the case for both the 'raw' 
values of aox etc. from the data and the 'intrinsic' values calcu- 
lated from de-reddened SED fits. This may imply that the optical 
threshold needs to be moved further into the far-UV before an Ed- 
dington ratio correlation is reported, if present. 

We also briefly discuss possible correlations of these spectral 
indices with mass. We remind the reader of the potential uncertain- 
ties with our K-band mass estimates discussed in S|3 .41 but we have 
found evidence that such estimates do trace more direct black hole 
mass estimators. We find only weak indications of a correlation 
between aband.X an d the black hole mass, but find that the corre- 
lation coefficient increases with decreasing reference wavelength, 
increasing from 0.30 (B-band) to 0.61 (UVW2 band; see Fig.[79t, 
if the outlier IRAS 09149-6206 (the AGN with the highest black 
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Figure 18. From left to right, top to bottom: Spectral indices connecting SED at the rest frame B-band, 2500A (ct x), U-band and UVW2-band with the SED 
at rest frame 2keV. Each spectral index is plotted against Eddington ratio. All conventions are as in Fig. 1171 



hole mass in this sample) is excluded. This relation implies that 
the spectrum flattens for higher mass objects, presumably indicat- 
ing that the fraction of UV emission directly from the accretion 
disc is less. A plot of bolometric correction against black hole mass 
(Fig. |20t displays this in a more pronounced fashion, confirming 
that the fraction of the luminosity appearing in X-rays, 1 / /12-iokeV 
(i.e. from the corona rather than the disc) increases with mass. This 
lends credibility to a scenario in which local, high mass black holes 
are radiating predominantly at low Eddington ratios: not only does 
the high mass lead to a low Eddington ratio, but lower bolometric 
correction reduces the bolometric luminosity further. Conversely, 
it is then expected that lower mass black holes exhibit higher Ed- 
dington ratios by virtue of their large bolometric corrections, in line 
with 'cosmic downsizing' scenarios in which the b ulk of accretion 
activi ty shifts to lower mass objects at low redshift dHeckman et al.l 
120041) . Since we are using the bulge luminosity to calculate the 
black hole mass, the possibility of a correlation of bolometric cor- 
rection with bulge luminosity is also interesting to consider. If real, 
it suggests that brighter bulges galaxies contain AGN accreting in 
states with a greater fraction of X-ray emission. This may indicate a 
link between the stellar populations of the host galaxy, the accretion 
state of the central black hole a nd the fuelling process; these the mes 
are discussed in more detail bv lKauffmann & Heckmaril d2008t) . 

The data from VF09 on the reverberation mapped sample do 
not exhibit any similar trends between either aox or bolometric 
correction and the black hole mass, despite the better quality mass 
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Figure 19. cko,b,u,UVW2-x against black hole mass. The key provides 
details of the symbols used. The best fit (solid line) excludes outlier IRAS 
09149-6206. 



estimates used. However, the reverberation sample is not well- 
selected and the quantities in VF09 do not include reddening and 
host galaxy corrections in the optical-UV, so do not necessarily 
provide a useful comparison for our purposes. 
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Figure 20. X-ray (2-10keV) bolometric correction against black hole 
mass. The best fit (excluding the outlier IRAS 09149-6206) has the form 
log(re 2 _iok c v) = 2 - 86 _ 0.1991og(M BH )/M Q (the spread in bolomet- 
ric corrections is a factor of ~ 2, by visual inspection). 



5.6 Host galaxy properties 

The host galaxy propert ies of the 9-month B AT catalogue are dis- 
cussed in some detail in IWinter et alj d2009h . Specific conclusions 
from their study include the lack of a convincing relation between 
host galaxy inclination and X-ray absorbing column Nu ; the large 
fraction of interacting systems in the whole catalogue (~54 per 
cent) and suggestions of evolution away from the elliptical, red 
hosts seen around AGN at z~l to systems with bluer colours. The 
accretion rate estimates in their paper do not suggest, however, that 
interactions/mergers are triggers for AGN activity, as the interact- 
ing systems do not exhibit a different distribution of accretion rates 
to the whole population. 

Our UVOT photometry results may provide interesting clues 
in understanding the host galaxy population of the catalogue. Here 
we present a discussion of the host-galaxy colours for our sample 
using the UVOT host galaxy magnitudes, and compare our findings 
with the literature (the more detailed work using KPNO data, Koss 
et al. in prep, will provide colours for the entire 9-month catalogue). 

There has been much previous work suggesting that AGN 
hosts occupy the so-called 'green valley' (SDSS g — r colour in 
the range ~ 0.5 — 0.9) in the colour-magnitude diagram, mid- way 
between the late-type 'red sequence' (g — r ~ 0.8 — 1.0) and star- 
forming 'blue cloud' (g — r ~ 0.2 — .5) galaxie s (see for example 
iNandra et alj2007l . lHickox et alj20"09T) . Recentlv. ls"chawinski et all 
J2009h present the colours for a small sample of 16 AGN from 
the BAT catalogue which have SDSS data and fin d that their hosts 
also li e in this intermediate regio n. The studies of lSilverman et al] 
J2008I) and lWestobv et al' 1 1200% are also of interest as they probe 
high and low redshift populations (from the Extended Chandra 
Deep Field South and the Sloan Digital Sky Survey, respectively). 
ISilverman et al] J2008h find that, accounting for biases caused by 
large-scale structure, AGN hosts at z > 0.8 have bluer hosts 
(U — V < 0.7) but the lo wer redshift counterpa rts (z < 0.6) tend 
to exhibit redder colours. IWestobv et al] d2007t> find a broad dis- 
tribution of colours for their composite AGN sample, towards the 
red end of the colour distribution for normal galaxies, and covering 
the green valley. We note t hat, in making our comparison with the 
results of I Silverman et al I < l2008l) . we do not have any biases with 
respect to large-scale structure due to the nature of the Swift-BAT 
catalogue, and additionally the only bias with respect to host galaxy 



properties would be the preferential selection of face-on galaxies 
(by virtue of our low Nh selection). 

Using the prescription in the UVOT CALDE0 for converting 
the UVOT magnitudes into standard Johnson magnitudes, along 
with the prescription for calculating SDSS colours from Johnson 
colours, we present our colour-magnitude diagrams in Fig.[2T] We 
present both the Johnson U — V colour against absolute V-band 
magnitude My and SDSS g — r colour against absolute SDSS r- 
band magnitude M(r); the g — r colour is calculated using the 
transformation g - r = 1.023(5 - V) + 0.16(£7 - B) - 0.187 
dKaraali et a"i]|2005l) and the r-band magnitude is interpolated us- 
ing the UVOT points. We provide U — V colours for 21 out of our 
26 objects for which both U and V host galaxy magnitudes were 
available from our GALFIT analysis; 3 of these objects did not 
have a B-band magnitude for the transformation to g — r, reducing 
our sample of objects to 18 for the latter plot. 

We see from Fig. [21] that the colours of our local sample 
sh ow a similar distributi on to the X-ray selected AGN in Fig. 3 
of lSilverman et al] J2008h . but w ithout a significant no ticeable red 
population. In comparison with IWestobv et al] J2007h . we find a 
distribution similar to their composite AGN sample, with a handful 
of very blue AGN. We ther efore seem to corroborate the sugges- 
tions in IWinter et al. I J2009h of an evolution away from red type 
host galaxies to bluer host galaxies, with the majority of objects 
in the 'green valley'. We also find som ewhat bluer colou r s than 
those found for the 16 BAT AGN hosts in lSchawinski et al] J2009I) 
although there are only five objects overlapping between the two 
samples (NGC 7469, NGC 5548, Mrk 766, MCG +04-22-042 and 
Mrk 1018; this list was obtained by private communication with the 
author). We caution that four of the bluest objects in Fig. 12 II (Ark 
120, Mrk 509, Mrk 841 and Mrk 352) have point-like morphologies 
in the UVOT images, and despite our efforts to minimize PSF mis- 
match, if such mismatch is present the wings of the PSF could be 
mistaken for the host galaxy giving rise to artificial 'hosts' which 
are influenced by the bluer colours of the AGN. These objects are 
depicted using crosses in Fig.|2T]for easy identification. This may 
indicate that the blue-ward shift is less pronounced than it may ini- 
tially seem, but the colours for the rest of our sample do lie in the 
blue end of the green valley. The more accurate and complete anal- 
ysis Koss et al. (in prep) will clarify these issues. 

We also caution that the PSF size of the UVOT may im- 
pose limitations on the accuracy of our host galaxy magnitudes. 
In tj3 - 1 - 1 1 we estimate the degree by which the nuclear flux is over- 
estimated with UVOT by comparison with HST ACS images for 
one source, NGC 4593. Based on the assumption that the flux spu- 
riously attributed to the nucleus should actually be attributed to 
the host galaxy, we can make a naive correction for this effect 
as follows. If the nuclear flux reported by UVOT is F^X° T = 
^iiuc T + ^spurious (for nuclear flux from HST F^f T and spuri- 
ous flux F S pu r ious), we can solve for -F spur i ous using the estimate 
of the error in the UVOT nuclear flux from CuTTl We can then 
obtain a better estimate of the host galaxy flux seen by HST, using 



F, 



HST 



galaxy 



= F. 



UVOT 



galaxy 



+ -Fspurious- The correction to the galaxy flux 



emerges as F spurious = F^ OT [l - (F^ u b c T /F^ OT )]. We em- 
ploy the fractions F^ U S C T /F^Z° T « 2.4/5.0 for the V-band and 
fa 1.2/1.8 for the U-band based on the analysis of NGC 4593, and 
calculate the corrections to the host galaxy magnitudes as follows: 



http://heasarc.nasa.gov/docs/heasarc/caldb/swift/docs/uvot/ 
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Figure 21. Left panel: Host galaxy colour U — V against absolute V-band magnitude for the host galaxies in the sample. Right panel: SDSS g — r colour 
against absolute r-band magnitude for the sample. Crosses represent objects with point-like morphologies, for which the host galaxy colours should be treated 
with caution. GALFIT returns errors of ± ~ 0.01 mag on magnitudes, but systematic errors due to the resolution of UVOT are likely to dominate (see Fig. 1221 
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It is not straightforward to predict whether the colours should 
generally increase or decrease under this correction, since the cor- 
rection is dependent on two things: the ratio of host galaxy flux to 
nuclear flux, and the degree of error in the nuclear flux intrinsic 
to UVOT, in a particular band. We plot the corrected U-V colours 
in Fig. [22] If the corrective factors assumed here are reasonable, 
we find that the U — V colours generally become bluer, in some 
cases substantially so. As expected, the corrected host galaxy V- 
band magnitudes are brighter due to the extra contribution from the 
nucleus. We do not make an attempt to correct the SDSS g — r 
magnitudes as these would require an estimate of the nuclear flux 
overestimate in the B-band as well, and the r-band magnitudes re- 
quire interpolation to wavelengths between the central UVOT filter 
wavelengths. It is therefore possible that better quality images may 
reveal even bluer hosts than our results indicate, but the limitations 
of the UVOT highlighted by this analysis suggest that the upcom- 
ing KPNO/SDSS work of Koss et al. (in prep) will provide more 
definitive estimates of the host galaxy colours. 



6 SUMMARY AND CONCLUSIONS 

We have presented SEDs for a well-selected subsample of local 
AGN from the Swift BAT 9-month catalogue of AGN. This study 
makes use of the unbiased and representative nature of the cat- 
alogue and offers numerous advancements over previous studies 
of the AGN SED. The simultaneous optical-to-X-ray observations 
from Swift used here provide a more accurate sn apshot of the ac- 
cretion energy budget in AGN, as pioneered by iBrocksopp et all 
and adopted by VF09. We are also able to increase confi- 
dence in measurements of the accretion luminosity by employing 
a subsample of AGN with minimal spectral complexity (no strong 
signs of partial covering or strong reflection) and lack of heavy 
absorption [log(A r H/cm _2 )<22] which could complicate the pic- 
ture: the X-ray spectra for these objects can all be fit with absorbed 
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Figure 22. Host galaxy colour U-V against absolute V-band magnitude, 
partially corrected for large PSF of UVOT. Key as in Fig. 12 II 



power-laws according to previous studies dWinter et al.ll2009h . al- 
lowing straightforward calculation of absorption-corrected X-ray 
luminosities. The data from Swift-UWOT typically cover at least 
four, if not all six available optical-UV filters on the instrument, 
and allow for a detailed reconstruction of the optical-UV SED. 
This is useful on two counts: firstly, in understanding the degree of 
host galaxy contamination in the observed SED for these objects; 
and secondly, in estimating how much dust extinction is responsible 
for the shape of the host-galaxy corrected nuclear SED. The inclu- 
sion of black hole mass estimates from K-band bulge luminosities 
allows accretion rates to be calculated using a uniform method for 
the sample. 

We consistently remove the host-galaxy contribution from the 
UVOT data points by fitting a 'PSF+disk' profile to the UVOT im- 
ages using GALFIT, and obtain a corrected nuclear flux in all avail- 
able filters. The host galaxy is often found to be a significant con- 
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taminant in the optical (V, B and U) bands, but far less significant 
in the UV bands. We also find that optical/UV emission lines are 
unlikely to contaminate the nuclear continuum significantly given 
the relatively large bandwidths of the UVOT filters (500-1000A). 
When fitted along with the simultaneous X-ray data with a sim- 
ple 'multicolour accretion disc + power-law' model, including both 
Galactic and intrinsic absorption, we are able to recover many im- 
portant parameters of the AGN SED, and produce an average model 
SED for lower-luminosity Seyferts with low accretion rates. The 
GALFIT fitting on the UVOT images also provides potential for 
understanding the host galaxy population of the unobscured BAT 
catalogue AGN. 

• In this representative sample, the 2-10 keV bolometric correc- 
tions cluster around values of 10-20. These bolometric corrections 
are essential for calculations of bolometric luminosities and accre- 
tion rates from X-ray observations, and therefore are of importance 
in the wider picture of scaling up the total X-ray energy density due 
to SMBH accretion in the X-ray background. 

• The Eddington ratios for the sample are typically below 0.1. 
Low bolometric corrections are expected for these low Eddington 
ratios (VF07, VF09), and furthermore, this may indicate that local, 
unobscured AGN are generally in an accretion state analogous to 
the iow/hard' state in Galactic black holes. This reinforces pre- 
vious findings that local Seyferts have lo w accretion rates, such 
as in early work bv lSun&Maikanl l ll989l) . who fit accretion disc 
models to optical-to-IR SEDs for a sample of 60 bright quasars 
and Seyferts. We arrive at the same conclusion using a different 
approach and employing a different sample: here we also use X- 
ray data in the calculation of accretion rates on a lower-luminosity 
sample selected for minimal absorption. 

• Some objects possess significant optical-UV dust redden- 
ing, even though these objects were selected for low absorb- 
in g gas column. This is broadly consist e nt wit h the observations 
of iKraemer et al.l ( l2000h . iBarcons et aTl d2003h and iMateos et all 
(2005) amongst others, and adds to the body of work suggesting 
that optical and X-r ay obscuration classifications can show sig nif- 
icant discrepancies llMaiolino et alJboOlL IComastri et alj|200l[ and 
ISilverman et alj|2005l are examples of studies which find evidence 
for the opposite effect in part or all of their samples; namely AGN 
with X-ray absorption but negligible optical-UV obscuration sig- 
natures). If we adopt a threshold of E(B — V) ~ 0.1 to divide 
heavily dust reddened and unreddened objects, ~ 5 — 6 out of our 
26 AGN (~ 20 per cent) are heavily dust reddened, consistent with 
the proportions of 'anomalous' AGN with mismatching optical and 
X-ray classifications found in some of the above studies, but this 
obviously depends on the threshold adopted. 

• The reddening-corrected ionizing luminosity fraction in- 
creases with Eddington ratio, as expected from combined accretion 
disc and corona models. 

• We confirm that spectral indicies linking the optical to the X- 
ray (qox etc.) show evidence for correlation with luminosity but 
not with Eddington ratio, corroborating the findings of previous 
studies. We find that extinction effects may need to be more care- 
fully considered when considering the correllation with luminosity, 
than has been done in previous works. De-reddening objects indi- 
vidually in constructing the aox — L • relation yields a weaker 

2500 A. 

correlation, which may need to be investigated in larger samples 
to check the robustness of this correlation and its use in constrain- 
ing the accretion process. The presence of a possible correlation of 
QUVW2X with mass, coupled with the attendant anti-correlation of 
bolometric correction with mass, may reinforce cosmic downsizing 



scenarios where lower mass black holes dominate AGN activity in 
the local universe by virtue of their larger accretion rates. 

• The host galaxies for this subsample generally lie in the green 
valley or blue-wards of it, in the standard colour-magnitude dia- 
gram. However, the UVOT images are of limited resolution and a 
better quality analysis with KPNO/SDSS (Koss et al. in prep) will 
shed more light on the host galaxy properties. 

• Adding the hard X-ray BAT data (14-195 keV) to these ob- 
servations can help with ascertaining the presence of a reflection 
hump in the spectrum. We selected four sources with the lowest 
hard X-ray variability, and find that three of these four showed a 
significant excess peaking at around ~30 keV which can plausibly 
be fit with a reflection hump. 



We have also discussed at length, the calibration of the black 
hole mass estimation method used to obtain these conclusions. The 
Eddington ratios and bolometric corrections could be significantly 
affected by these uncertainties. We estimate the degree of this ef- 
fect by comparison with results using reverberation mass estimates. 
This analysis reveals somewhat higher bolometric corrections and 
Eddington ratios, but the bulk of the sample lies consistently at 
lower Eddington ratios and exhibits low bolometric corrections. 
The bulge luminosity-based and reverberation mapping methods 
require more accurate calibration before more definitive statements 
on the magnitudes of these uncertainties can be constrained. 

The hard X-ray selection criteria used to specify the BAT cat- 
alogue eliminate potential biases. As an illustration of this, we no- 
tice the presence of substantially more low accretion rate objects in 
this sample than the samples of VF07 and VF09 (the former was 
UV selected by FUSE, the latter selected for being optically bright 
enough for reverberation mapping). The low accretion ra tes found 
here broadly corroborate the results o flWinteretafl ( l2009h using X- 
ray luminosities (with no bolometric correction applied). Although 
these objects have low-absorption and are situated at low redshift, 
they are similar in power to the obscured, higher redshift sources 
thought to be r esponsible for the bulk of the X-ray background. 
As suggested in iFabiarj 120041) , we confirm that a lower bolomet- 
ric correction is appropriate for these objects, based on our multi- 
pronged approach to recovering the true bolometric luminosities in 
these sources. Bolometric corrections of the magnitude seen here 
are appropriate for reconciling the energy density from the X-ray 
background with the local black hole density. It remains to be seen 
if such properties are displayed in large samples of obscured AGN, 
the dominant class of AGN contributing to the X-ray background. 

Useful extensions to this work would be to acquire simul- 
taneous Swift observations for the remaining objects in the BAT 
catalogue, which would potentially confirm the trends we identify 
here, or highlight unknown biases in our chosen sample. What is 
of particular importance next is to understand the workings of the 
central engine of obscured AGN. In a companion paper (Vasude- 
van et al. in prep) we hope to address this issue for both obscured 
and unobscured AGN in the Swift/BAT catalogue, using the re- 
processed IR emission along with the hard X-ray BAT observa- 
tions to estimate their bolom etric properties, analogous to the work 
done bv |Pozzietai] ( l2007l) for eight higher redshift quasars. The 
well-studied 9-month catalogue again provides an excellent start- 
ing point for such work, and follow-up studies on the 22-month 
catalogue dTueller et al]|2009h will offer new scope for constrain- 
ing the accretion physics of a larger, representative sample of AGN. 
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Table Al. Results from SED fitting. 

(1) - Log of central BH mass from 2MASS K-band magnitude. Errors from K-band magnitudes translate into errors of ±0.1 in values for log(MBH/MQ), 
but we refer the reader to the systematics discussed in the text. 

(2) - Photon index from fit to 0.3-10keV regime. 

(3) - Log of 2-10keV luminosity from power-law fit to 0.3-10keV regime. 

(4) - Spectral index q x calculated by interpolation between available UVOT data points to determine the 2500A luminosity. 

(5) - Intrinsic extinction E(B - V) from fitting ZDUST(DISKPN) model combination to UVOT data. 

(6) - Log of bolometric (0.001-100keV) luminosity, corrected for X-ray absorption and optical-UV dust reddening. 

(7) - Spectral index q x calculated from the full optical-to-X-ray model fit, corrected for optical-UV dust reddening. 

(8) - Eddington ratio. 

(9) - Hard X-ray 2-10keV bolometric correction £bol/^2-l0keV- 

Random error estimates are provided on values for intrinsic extinction, Eddington ratio and bolometric correction; random errors on all other quantities are 
omitted as systematic errors generally dominate. 
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